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ABSTRACT 

We investigate the evolution of the faint-end slope of the luminosity function, a, using semi- 
analytical modeling of galaxy formation. In agreement with observations, we find that the slope 
can be fitted well by a{z) = a + bz, with a = —1.13 and b = —0.1. The main driver for the evolution 
in a is the evolution in the underlying dark matter mass function. Sub-L* galaxies reside in dark 
matter halos that occupy a different part of the mass function. At high redshifts, this part of the mass 
function is steeper than at low redshifts, and hence a is steeper. Supernova feedback in general causes 
the same relative flattening with respect to the dark matter mass function. The faint-end slope at 
low redshifts is dominated by field galaxies and at high redshifts by cluster galaxies. The evolution of 
a{z) in each of these environments is different, with field galaxies having a slope b = —0.14 and cluster 
galaxies b — —0.05. The transition from cluster-dominated to field-dominated faint-end slope occurs 
roughly at a redshift ~ 2, and suggests that a single linear fit to the overall evolution of a{z) might 
not be appropriate. Furthermore, this result indicates that tidal disruption of dwarf galaxies in clus- 
ters cannot play a significant role in explaining the evolution of a{z) at z < z*. In addition we find 
that different star-formation efficiencies a* in the Schniidt-Kennicutt-law and supernovae-feedback 
efficiencies e generally do not strongly influence the evolution of a{z). 

Subject headings: galaxies: evolution ~ methods: numerical 



1. INTRODUCTION 

The galaxy luminosity function (LF) is one of the 
cornerstones in our understanding of galaxy formation 
and evolution. Since the introduction of a fitting func- 
tion for its shape by ISchechteil ()1976f ) , the origin of the 
form of the LF function has been a powerful constraint 
onmodel building (e.g.. iBenson et aDl2003l : ISamui et al.l 
l2007f ). While recent work has focused som ewhat on the 
luminous end, its evolution with redshift (B rown et al. 
2007.) and the role of dr y mergers (.Khochfar fc Burkert 
200a iNaab et all l2006( l. the faint-end provides addi- 



tional important clues on galaxy formation. Systematic 
studies of the faint-end slope in the local universe re- 
veal dif ferences be tween high and low density environ- 
ments (jTrenthaml [T998il . as well as for galaxy samples 
split by morphologies (e.g. iMarzke et al.l I1994D . The 
underlying physical processes that shape the faint-end 
of the LF are generally associated with feedback from 
supernovae that is effective in heating ga s and driving 
wind s in shallow gravitational potentials (|Dekel fc Silkl 
Il986l ). Although the implementation of supernova feed- 
back in galaxy f ormation models ha s been extensively in- 
vestigated fe.g. iBenson et al]|2003f) for the local galaxy 
LF, its impact on the redshift evolution on the faint- 
end has not been as well studied. The recent wealth 
of LF measured t o very faint niagnitudes in the rest- 
frame B-band (e .g. iBlanton et al."2003'; ' Wolf et al.ir2003t 
iMarchesini et al . 2007; Rvan ct al. 2007) and the rest- 
frame FUV (e.g.lY an & Windhorsti 120045 IWvder et aP 
l2005t iBouwens et al. 2006) allows us to test models with 
high accuracy. 

The purpose of this letter is to investigate the under- 
lying driving mechanism for the redshift evolution of the 
faint-end slope. Furthermore, we investigate the impact 



of supernova feedback on the rate of star-formation by 
varying the relevant efficiency parameters. 

2. MODEL 

In the following we briefly outline our basic mod- 
elling approach and refer the reader for more details to 
Eiochfar & Burkert (2005), Khochfar & Silk (2006) and 
reference therein. We generate merger trees for dark mat- 
ter halos using a Monte-Carlo appr oach based on the ex- 
tende d Press- Schechter formalism (jSomerville fc KolattI 
[19991) . As we aim to model the faint-end of the LF 
to high redshifts, we need to make sure that the mass 
resolution in our simulations is sufficient. We generate 
merger trees from dark matter mass functions between 
z = and z = 6, and find that resolving each individual 
merger tree down to a mass resolution of M^^^^^ = 5x10^ 
Mq and Mj^ijj = 10^ Mq at z < 3 and z > 4, re- 
spectively, gives robust results. Once a tree reaches 
Mmin, we start moving the tree forward in time includ- 
ing physical processes associated with the baryons within 
each dark matter halo that include gas cooling, star- 
formation, supernova feedback, reionization and merg- 
ing of galaxies on a dynamical friction time-scale. As 
the focus of this letter is on the faint-end of the lumi- 
nosity function, we will omit includin g prescriptions for 
AGN-feedback (e.g. 'Bower et al."2006) or environmental 
effects (Khochfar & Ostriker 2007) that mainly influence 
the bright-end of the luminosity function. 

The largest impact on the slope at the faint-end 
comes fro m star- formation and associated supernova 
feedback (|Dekel fc Silkl [198^ . Faint galaxies generally 
occupy small dark matter halos with shallow potential 
wells, that allow effective reheating of cold gas in the 
ISM by feedback from supernovae. We model star- 
formation in galaxies using a parametrisation of the 
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global Schmidt-Kennicutt-law (jKennicuttl 1 1 998f ) accord- 
ing to which = a*Mcoid/tdym where a* is a free 
parameter that is indicative of the efficiency of star- 
formation, Mcoid is the mass in cold gas in the galac- 
tic disk, and tdyn is the dynami cal time-scale of th e 
galaxy. Following the arguments by lDekel fc Silk! (|1986f ). 
we model the amount of cold gas reheated by feedback 
from supernovae with Msn = 4:eMt:r]sNEsN /SVj^g^^, 
with e as a free parameter that controls the feedback 
efficiency, rjsN the number of supernovae per solar mass 
of stars formed, Esn = 10^^ erg is the energy released 
by each supernova, and Vmax is the maximum circular 
velocity of the dark matter halo in which the galaxy re- 
sides. 

For each individual galaxy in our simulation, we store 
the star-formation history and generate its i3— band and 
FUV rest-frame luminosity function at various redshifts 
using the stellar population models of lBruzual fc CharlotI 
([2b03,) . The faint-end of the luminosity function is then 
fitted by a simp le power law with slope a as defined in 
ISchechteij()1976[ ). We fit the faint-end LF at each redshift 
with a power law spanning a range of four magnitudes 
at z < 3 and at least two magnitudes at 2; > 3, start- 
ing at the lowest magnitude Lmin that is unaffected by 
the mass resolution of the simulation. We choose this 
approach over fitting the whole LF with a Schechter-fit, 
because we are missing physical effects (see below) in 
our model that are responsible for shaping the bright- 
end and the knee of the LF. In addition, we increase 
the number of magnitude bins and make sure that the 
fitted values for a are unaffected by the bin-size. Fig.l 
shows the luminosity function at magnitudes larger than 
the corresponding minimum magnitude for Lmm- In this 
study, we simulate a volume of 10® Mpc^, which allows 
us to calculate a robustly up to a redshift z ~ 6. 

Throughout this paper, we use the following set of cos- 
mological parameters based on the three year Wilkinso n 
Microwave Anisotropy Probe data (jSpergel et al.ll2007l ): 
flo = 0.27, = 0.73, nt/no = 0.17, erg = 0.77 and 
h = 0.71. 

3. RESULTS 

There is significant evidence that the faint-end slope of 
the galaxy luminosity function shows a measurable de- 
pendence on redshift, which can be fitted by a linear law 
of the form a(z) = a + bz, with a between —1.12 and 
— 1.17, and b between —0.12 and —0-11 (for recent obser- 
vations. seelSawicki fc Thompson|[2006tlMarchesini et al.l 
120071: iRvan et al.ll2007D . WithhT the hierarchical struc- 
ture formation paradigm, one naturally expects such be- 
haviour, considering that the slope of the dark matter 
mass function below Mdm,* is aoM ^ 2 and that the 
objects that form in these halos continue to grow by 
continued star-formation a nd mergers with each other 
(jKhochfar fc Burkertll200ll ). hence flattening the slope. 
Fig. 1 shows the predicted model luminosity function 
at various redshifts in the rest-frame B and UV. The 
simulated and observed luminosity functions are in fair 
agreement at the faint-end. The luminous end however, 
shows deviations at low redshift which are due to missing 
feedback sources in massive galaxies such as AGN or to 
environmental effects. In Fig. 2, we show the predicted 
evolution of a{z) for our best-fit local model. The free 
parameters a* and e in this model are chosen to give the 
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Fig. 1. — Comparison between model LFs with e = 0.6 and 
a = 0.02 (symbols), and Schechter fits to various observations of 
the LF at different redshifts (solid lines). The dashed lines show 
our power law fits to the faint-end of the LF. Top left panel: Rest- 
frame B— band LF from Norbcrg et al. (2002), where we assumed 
that B = bj+0.12 mag. The shaded area indicates the region of the 
LF that is not well matched due to missing feedback effects. Top 
right panel: Rest-frame _B— band LF from (Marchesini et al. 2007|). 
Bottom left: UV-LF fr om Bouwens et al. (2007). Botto m right: 
UV-L F from lYan &: W indhorst (2004b). Following Bouw ens et al.l 
(120061) . we assume an average dust correction 0.4 mag at z > 4. 

best fi t to various local observations fsee lKhochfar &: Silkl 
l2006f ). For consistency with the majority of observations, 
we calculate the faint-end slope for the rest-frame FUV 
at z > 4 and at lower redshifts for the rest-frame i?-band. 
We indeed find an evolution in a with redshift that is in 
fair agreement with the observed evolution. 

The immediate question that arises is, what influences 
and is the main driver for the evolution in a? Generally, 
supernova feedback is considered the dominant mecha- 
nism in shaping the faint-end of the lum inosity function 
(|Dekel fc Silk|[l986l : llmson et al.l[200l . The shaded re- 
gion in Fig. 2 shows the range of linear fits to a{z) that 
we find by varying the star-formation efficiency between 
a* = 0.02—0.1 and the supernovae feedback efficiency be- 
tween £ = 0.2 — 0.6. We infer only a very modest change 
in a(z) for reasonable choices of feedback efficiencies, and 
therefore conclude that another process must be respon- 
sible for the observed evolution in a{z). 

The mass function of dark matter halos is known 
to show a strong e volution with redshift (e.g. 
iPress fc Schechteil 11974) . The galaxies contributing to 
the luminosity function around are mostly central 
galaxies in their dark matter halo s, i.e. the most lumi- 
nous galaxy within the halo (e.g. iKhochfar fc Ostriked 
I2OO7I ). It is therefore not unreasonable to assume a con- 
nection between the evolution of a{z) and that of the 
dark matter mass function. When considering the lumi- 
nosity of central galaxies residing in dark matter halos 
of the same mass at different redshifts, we find that at 
early times, central galaxies are up to three magnitudes 
brighter than their counterparts in low redshift halos. 
This is even the case for halos hosting sub-L,^ galaxies. 
Similar results have been reported by iKobavashi et all 
|2007i) . who showed that dwarf galaxies at early times 
are not affected by supernova feedback in their Simula- 
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Fig. 2. — The slope a at different redshifts as predicted by the 
best fit local model with a* = 0.02, e = 0.6. Filled symbols show 
results from the simulation and the solid line is the best fit to the 
simulation data. Errorbars indica te l-a errors . The dashed lines 
show the fit to the compiled data in lRvan et al.l ((2007). The shaded 
region shows the range of linear fits to ce{z) that we find when 
varying the star-formation and supernovae feedback effici encies as 
discussed in the text. Stars are the compilation from Ryan et al.l 
II2007IV The filled squares and triangles show recent results from 
lOesch et al.l ([20071) and Bouwcns et al. (2007), respectively. 

tions because cooling times are very short in these halos 
and the energy injected by the supernovae is rapidly dis- 
sipated away. The slope in the region of dark matter ha- 
los that host sub-L* galaxies is steeper at high redshift, 
and consequently so is a. The same is true for other 
choices of a, and e, thereby explaining why we do not 
find any strong dependence of a{z) on these parameters. 
It should be noted however, that modelling these param- 
eters with a strong redshift dependency will be able to 
enhance or weaken the evolution of a 

We continue analyzing the evolution in a{z) by dis- 
tinguishing between cluster and field galaxies and their 
relaxation to the overall luminosity function at the faint- 
end. In Fig. 3, we present a{z) for progenitor galaxies 
of present-day cluster and field galaxies from our simu- 
lations. Here we define cluster environments by present- 
day dark matter halos above 10"'^'* Mq, and field envi- 
ronments by halos with masses below 10^^ M©. In field 
environments a is steeper and evolves more strongly than 
in cluster environments. At early times, the first galaxies 
to appear are most likely in high— cr fiuctuations, which 
will result in present-day galaxy clusters. Consequently, 
the faint-end luminosity function at high redshifts will 
be dominated by present-day cluster members, and the 
faint-end-slope of the overall galaxy population at high 
redshifts is flatter than that for the field luminosity func- 
tion alone at the same redshift. When considering the 
relative weight of field galaxies to the overall galaxy 
population at the faint-end, one can estimate the red- 
shift at which the transition from cluster-driven to field- 
driven evolution in a occurs. We find that this transition 
roughly occurs at ~ 2 with a slight dependency on the 
definition of environment. 

4. DISCUSSION AND CONCLUSIONS 
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Fig. 3. — The slope a at different redshifts as predicted by the 
best fit local model. Filled circles show results for typcal field 
environments and filled stars show the results for cluster environ- 
ments, as defined by their present day dark halo mass. The solid 
and dashed line are the fits to the modeled evolution in the field 
and cluster, respectively. Here we define cluster environments by 
present-day dark matter halos above 10^* Mq and field environ- 
ments, by halos with masses below 10^2 Mq. 

In this Letter, we presented predictions for the redshift 
evolution of the faint-end slope of the luminosity func- 
tion within the ACDM-scenario. In general, we find the 
same trend as in the recent observations i.e., a steepening 
faint-end slope a with redshift, which can be well-fitted 
by a simple linear fit a{z) = a + bz where the observa - 
tions find a ~ 1.17 and h ~ -0.11 (iRvan et al.ll2007| ). 
Our simulations predict a ~ 1.13 and b ~ —0.1, in good 
agreement with the observations considering the large 
uncertainties, especially at high redshifts. 

Our simulations confirm previous results that the fiat- 
tening of the faint-end slope a with respect to the slope 
in the dark matter mass function can be well explained 
by supernova feedback. However we additionally show 
that a is steeper at higher redshifts mainly, due to the 
dark matter mass function being steeper for the range of 
halo masses that host sub-L* galaxies, suggesting that 
the evolution of a traces closely that of the underlying 
dark matter mass function. 

The contribution of the progenitor population of 
present-day field and cluster galaxies plays a significant 
role in shaping the evolution of a . For field galaxies, 
the evolution of a{z) is stronger, with b = —0.14 and 
a = —1.16 than for cluster galaxies. In our simulations, 
we find that at redshifts z > 2, the faint-end is domi- 
nated by galaxies ending up in present-day clusters. This 
transition redshift is dependent on the value of dg, which 
normalizes the power spectrum and regulates the redshift 
at which structures of a given mass typically form. Ad- 
ditionally, the slope of the fiuctuation spectrum at small 
scales will influence z*. Precise high redshift measure- 
ments of the contributions of these two populations to 
the faint-end of the luminosity function in future surveys 
with e.g. the HST WFC3 will help to pin down z*. One 
potential problem for future surveys will be a possible 
bias towards cluster gal axies as they might e xperience 
induced star- formation (|Marcillac et al.ll2007f ). increas- 
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ing their surface brightness and making them more eas- 
ily detectable. This effect will shift ^* to higher redshifts 
and needs to be taken into account carefully. Obser- 
vational selection effects will affect the observed faint- 
end LF-slope in Fig. 2. Some observational selection 
effec ts (i.e. catalogue incoin pleteness and natural confu- 
sion, iWindhorst etani200l can make the observed faint- 
end slope flatter than the true one, while others (e.g., 
SB-dimming) could make the observed slope somewhat 
steeper than the true one, depending on the exact in- 
trinsic object size distribution. A number of groups cor- 
rect for incompleteness either through MC-simulations 
(e.g. lYan fc Wi ndhorst* '2004b') or through cloning tech- 
niques (e.g. iBo uwcns ct al. 2006), and find similar faint- 
end slopes when following different procedures. When 
judging the data, however, one must keep these observa- 
tional biases in mind. Ultimately, these issues can only 
be resolved with deeper JWST data to AB=31-32 mag. 

Tidal disruption of dwarf g alaxies in clusters as seen in 
high resolution simulations (iTormen et al.|[l"998f ) can in 
principle change the slope a. Our results, however, sug- 
gest that at a transition redshift of = 2, the evolution 
of a changes from being dominated by cluster galaxies 
to being dominated by lield galaxies. It is therefore not 
likely that a large amount of the evolution in a at z < 
is driven by tidal disruption of faint galaxies. An addi- 
tional implication from the transition at is that the 
evolution of a{z) is better fit by a linear function with a 
break at z,. 

The flattening of the slope a with respect to the un- 
derlying dark matter slope suggests the interesting pos- 
sibility of estimating the time-scale over which super- 
novae operate. Assuming that the firs t Pop II stars were 
formed sometime before reionization (|Yan &: Windhorst 
I2004U) . and that supernovae type la originate roughly 
< 1 — 2 Gyr after the bulk of the first Pop II star- 
formation, one would expect an increase in energy in- 



jection into the interstellar medium at a redshift corre- 
sponding to this time lag. This additional energy input 
will hinder star-formation and contribute to a further in- 
crease in the mass-to-light ratio of galaxies and hence 
to an even stronger flattening of the slope a. It will be 
crucial to have robust measurements of a over a wide 
range of redshifts to probe the onset of the first signifi- 
cant feedback contribution from type la supernovae. Fur- 
thermore, probing the faint-end slope at redshifts z > 6, 
before the significant onset of type II supernovae will al- 
low us to measure the underlying dark matter slope very 
accurately. 

Our approach has certain shortcomings. The model 
presented here did not include any time delay prescrip- 
tions for the various SN types, but instead assumed 
instantaneous feedback. More detailed modeling of the 
time delays and its influence on the faint-end slope 
will be presented elsewhere (in preparation). Our 
treatment of supernovae feedback is very simplistic, and 
more detailed hydrodynamical simulations including 
a multi-phase medium will show if this general trend 
which we report can be recovered. First generations 
of such simulations indeed show that SN type II 
that are generated in dense star clusters explode into 
bubbles of hot gas and are therefore less efficient at 
feedback into the ISM (Mac Low & Ferrara 1999|) 
compared to SN type la, which go off at random places 
in the galaxy and, can have more effect on the early ISM. 

We would like to thank Scth Cohen, Evan Scanna- 
pieco and Richard Bouwens for helpful discussions and 
the anonymous referee for his useful comments. This 
work was supported by HST grants HST-GO-10530.07 
(to RAW) and HST- AR- 109 74. 01 (to RER) from STScI, 
which is operated by AURA for NASA under contract 
NAS 5-26555, and by NASA JWST grant NAG 5-12460 
(to RAW). 



REFERENCES 



Benson, A. J., Bower, R. G., Frenk, C. S., Lacey, C. G., Baugh, 

C. M., & Cole, S. 2003, ApJ, 599, 38 
Blanton, M. R., et al. 2003, ApJ, 592, 819 

Bouwens, R. J., lUingworth, G. D., Blakeslee, J. P., & Franx, M. 

2006, ApJ, 653, 53 
Bouwens, R. J., lUingworth, G. D., Franx, M., & Ford, H. 2007, 

ArXiv e-prints, 707, arXiv:0707.2080 
Bower, R. G., Benson, A. J., Malbon, R., Helly, J. C, Frenk, C. S., 

Baugh, C. M., Cole, S., & Lacey, C. G. 2006, MNRAS, 370, 645 
Brown, M. J. I., Dey, A., Jannuzi, B. T., Brand, K., Benson, A. J., 

Brodwin, M., Croton, D. J., & Eisenhardt, P. R. 2007, ApJ, 654, 

858 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Cole, S., et al. 2001, MNRAS, 326, 255 
Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

Huang, J.-S., Glazebrook, K., Cowie, L. L., & Tinney, C. 2003, 
ApJ, 584, 203 

Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, 
S., Evrard, A. E., Couchman, H. M. P., & Yoshida2007, N. 2001, 
MNRAS, 321, 372 

Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 

Khochfar, S., & Burkert, A. 2001, ApJ, 561, 517 

Khochfar, S., & Burkert, A. 2003, ApJ, 597, L117 

Khochfar, S., & Burkert, A. 2005, MNRAS, 359, 1379 

Khochfar, S., & Silk, J. 2006, MNRAS, 370, 902 

Khochfar, S., & Ostriker, J. P. 2007, ArXiv e-prints, 704, 
arXiv:0704 .2418l 

Kobayashi, C, Springel, V., & White, S. D. M. 2007, MNRAS, 
376, 1465 

Marchesini, D., et al. 2007, ApJ, 656, 42 

Marcillac, D., Rigby J. R., Rieke, G. H., & Kelly D. M. 2007, ApJ, 
654, 825 



Marzke, R. O., Geller, M. J., Huchra, J. P., & Corwin, H. G., Jr. 

1994, AJ, 108, 437 
Mac Low, M.-M., & Ferrara, A. 1999, ApJ, 513, 142 
Naab, T., Khochfar, S., & Burkert, A. 2006, ApJ, 636, L81 
Norberg, P., et al. 2002, MNRAS, 336, 907 

Oesch, P. A., et al. 2007, ArXiv e-prints, 706. larXiv:0706.2653l 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Ryan, R. E., Jr., et al. 2007, ArXiv Astrophysics e-prints, 

arXiv:astro-ph/0703743 
Schechter, P. 1976, ApJ, 203, 297 
Sawicki, M., & Thompson, D. 2006, ApJ, 648, 299 
Samui, S., Srianand, R., & Subramanian, K. 2007, MNRAS, 377, 

285 

Schiminovich, D., et al. 2005, ApJ, 619, L47 

Somerville, R. S., & Kolatt, T. S. 1999, MNRAS, 305, 1 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Tormen, G., Diaferio, A., & Syer, D. 1998, MNRAS, 299, 728 

Trentham, N. 1998, MNRAS, 294, 193 

Trentham, N., Sampson, L., & Banerji, M. 2005, MNRAS, 357, 783 
Wolf, C, Wisotzki, L., Borch, A., Dye, S., Kleinheinrich, M., & 

Meisenheimer, K. 2003, A&A, 408, 499 
Windhorst, R. A., Hathi, N. P., Cohen, S. H., Jansen, R. A., 
Kawata, D., Driver, S. P., & Gibson, B. 2 007, Advances in Sp ace 
Research, Vol. 9163, p.001-013, in press l |astro-ph/0703171| i 
Wyder, T. K., et al. 2005, ApJ, 619, L15 
Yan, H., & Windhorst, R. A. 2004a, ApJ, 612, L93 
Yan, H., & Windhorst, R. A. 2004b, ApJ, 600, LI 



